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Abstract 

The reactions pp — > pf(X°)p s , where X° is observed decaying to 7r 7r 7r 7r° , 7r + 7r~7r + 7r~ and 
7r + 7r _ 7T 7r° , have been studied at 450 GeV/c. There is evidence for an ci2(1320)7r decay mode 
of the ?72(1645) and 772 (1870) in the 7r + 7r~7r°7r° and 7t + 7t~7t + tt~ final states. The ^(1950) is 
consistent with being a single resonance with a dominant /2(1270)7T7t decay mode. The /o(1370) 
is found to decay dominantly to pp while the /o(1500) is found to decay to pp and ao. 
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The WA76, WA91 and WA102 collaborations have studied the centrally produced 7t + 7t~7t + tt~ 
final state in the reaction 

pp — > pf(n + n^7Y + n^)p s (1) 

at 85 300 0, and 450 GeV/c []3], |]. The subscripts / and s indicate the fastest and slowest 
particles in the laboratory respectively. In addition to the /i(1285), which was observed at all 
energies, peaks were observed at 1.45 and 1.9 GeV in the 300 and 450 GeV/c data. In contrast, 
no clear evidence was seen for these states in the 85 GeV/c data of the same experiment 0. 
The increased prominence of these states with increased incident energy is consistent with 
the formation of these states via a double Pomeron exchange mechanism, which is predicted to 
be a source of gluonic states p. 

The peak at 1.9 GeV, called the / 2 (1950), was found to have I G J PC = 0+2++ and decay 
to /2(1270)7T7r and <22(1320)7r 0, f3|. It was not possible to determine whether this was one 
resonance with two decay modes, or two resonances. However, in a recent analysis of the 
centrally produced r]7i + 7i^ final state by the WA102 experiment ||, no evidence was found for 
a J PC = 2 ++ a 2 (1320)7r decay mode. The peak at 1.45 GeV was shown to have I G J PC = 0+0++ 
and decay to pp. It has been described as being due to the interference between the /q(1370) 
and the /o(1500) || £§] which implies that both states should have a substantial pp decay mode. 
However, the initial analysis of the 4n channel in pp by Crystal Barrel showed that the 
/o(1500) decays via oo with a very weak coupling to pp. A recent preliminary analysis |§ of 
the Crystal Barrel data concludes that the /o(1370) decays strongly to oo and pp while the 
/o(1500) decays to oo and 7i*tt where the 7r*7r decays to pir . 

Therefore, there are still several unanswered questions on the decays of the /2(1950), 
/o(1370) and /o(1500). In order to address these questions, in this paper an analysis of several 
different 4n decay modes will be presented. In particular, to investigate the oo contribu- 
tion in the decays of the / (1370) and / (1500). One of the major problems of studying the 
7r + 7T~7r + 7r~ final state is the number of possible isobar decay modes that are present. Therefore 
a study will be made firstly of the 7r 7r 7r 7r final state which has the nice feature that it is 
free from any contribution from p(770), ai(1260) or a2(1320) isobars. Next an analysis will be 
presented of the 7r + 7r~7r 7r° decay mode which has different systematic effects compared to the 
7r + 7r~7r + 7r~ decay. Finally a reanalysis of the 7r + 7r~7r + 7r~ decay mode will be presented. 

The data come from the WA102 experiment which has been performed using the CERN 
Omega Spectrometer, the layout of which is described in ref. ||. The reaction 

P ^ P/ (AW) Ps (2) 

has been studied at 450 GeV/c. Reaction (||) has been isolated from the sample of events having 
two outgoing charged tracks and eight 7s reconstructed in the GAMS-4000 calorimeter, by first 
imposing the following cuts on the components of missing momentum: | missing P x \ < 14.0 
GeV/c, (missing P y \ < 0.20 GeV/c and (missing P z \ < 0.16 GeV/c, where the x axis is along 
the beam direction. A correlation between pulse-height and momentum obtained from a system 
of scintillation counters was used to ensure that the slow particle was a proton. 

Fig. [Tja) shows the two photon mass spectrum for 87-events when the mass of the other 
three 27-pairs lies within a band around the 7r° mass (100-170 MeV). A clear 7r° signal is 
observed with a small background. Events belonging to reaction (H) have been selected using 
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a kinematical fit (8C fit, four-momentum conservation being used and the masses of four 7r°s 
being fixed). Events corresponding to the reaction i]7r°, where the r) decays to 3tt°, have been 
removed by requiring M(3n°) > 0.6 GeV. The major background to the 7r°7r 7r 7r° final state 
comes from the decay of the rf and /i(1285) to i]tt°tt° where the t] decays to 37r° and one ir° 
is undetected. This produces events in the threshold to 1.3 GeV region of the 7r 7r 7r 7r mass 
spectrum. A kinematical fit has been used to remove these events. Using real rjirir events a 
simulation has been performed to estimate the residual contamination in the 7r°7r 7r 7r° final 
state which is found to be less than 5 %. 

The resulting 7r 7r 7r 7r° effective mass spectrum is shown in fig. |l|b) and consists of 1438 
events. There is a broad enhancement in the 2 GeV region, compatible with the /2(1950), but 
there is no clear peak in the 1.45 GeV region. 

One of the aims of this paper is to investigate the role of the a in the decays of the 4ir 
final state. In order to determine the amount of aa a parameterisation of the tttt S-wave (a) is 
required. Parameters for the a may be process dependent. The nature of the a is not clearly 
established and parameterisations of it may be process dependent. It may be a resonant state 
or just two pions produced in relative S-wave. However, the parameterisation of the tttt S-wave 
is crucial in all these analyses in order to determine the amount of aa observed. Assuming 
that the 1.3—1.5 GeV 7r 7r 7r 7r mass region is dominated by the decay to aa, then in order to 
determine the best parameters for the a, a study has been made of the 7r°7r° mass spectrum for 
the region 1.3 < M(7r 7r 7r 7r° ) < 1.5 GeV, shown in fig. |l]c). A Monte Carlo simulation has 
been performed using the decay X° — > aa with a — > 7r°7r°. Formally, we describe the shape of 
the a using a spin relativistic Breit-Wigner amplitude of the form: 



In order to describe the centrally produced mass spectrum the resulting function ^(M^)! 2 
has been multiplied by the kinematical factor (M„ — 4m 2 r ) 1//2 /M| 7r flnj. The mass and width 
of the a have been varied and the Monte Carlo data passed through the detector simulation 
and the resulting mass distributions were compared to the real data. The parameters that best 
describe the data are m = 450 ± 50 MeV and T = 600 ± 50 MeV and the resulting mass 
spectrum (solid histogram) is superimposed on the real 7r°7r° mass spectrum shown in fig. |]c). 

A spin-parity analysis of the 7r°7r 7r 7r° channel has been performed in 120 MeV 7r 7r 7r 7r° mass 
bins using an isobar model and the method described in ref. [[J, assuming that only the aa 
and /2(1270)7T7r intermediate states with J PC = ++ or 2 ++ contribute. The data can be fitted 
using only the J PC = ++ a a wave and the J PC = 2 ++ /2(1270)7T7t wave with Jz = which 
are shown in figs. |l|d) and e) respectively. Superimposed on fig. |l]d) is the distribution that 
would be expected for the /o(1500); as can be seen there is little evidence for any /o(1370) 
contribution. Superimposed on fig. |]e) is what would be expected from the /2(1950); as can be 
seen it well describes the J PC = 2 ++ /2(1270)7T7t wave, showing that the peak at 1.95 GeV is 
consistent with being a single resonance. 

An analysis has next been performed on the centrally produced 7r + 7r~7r 7r° system. The 
reaction 




pp — > pf(n + n 7T 7r°)p, 



's 



(3) 
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has been isolated from the sample of events having four outgoing charged tracks and four 
7s reconstructed in the GAMS-4000 calorimeter, by first imposing the following cuts on the 
components of the missing momentum: |missing P x \ < 17.0 GeV/c, |missing P y \ < 0.16 GeV/c 
and |missing P z \ < 0.12 GeV/c. The two photon mass spectrum (not shown) when the mass of 
the other 27-pair lies within a band around the 7r° mass (100-170 MeV) shows a clear 7r° signal 
with small background. Events belonging to reaction (^) have been selected using a kinematical 
fit (6C fit, four-momentum conservation being used and the masses of two tt°s being fixed). 
Events containing a fast A ++ (1232) were removed if M(pj7r + ) < 1.3 GeV, which left 283 408 
centrally produced 7r + 7r~7r°7r° events. 

Fig. shows the 7r + 7r~7r°7r° effective mass spectrum. The mass spectrum is very similar 
to that observed previously in the 7r + 7r~7r + 7r~ channel Q, namely, a clear peak at 1.28 GeV 
associated with the /i(1285), a peak at 1.45 GeV and a broad enhancement around 2.0 GeV. 

A spin-parity analysis of the 7r + 7r~7r°7r° channel has been performed using an isobar model 
and the method described in ref. |3j. Assuming that only angular momenta up to 2 contribute, 
the intermediate states considered are 

aa, up , 

P + P~, P ± (^^°)swave, P ± (^71°) Pwave , (n T 7Y°) Dwave , 

ai(1260)7r, a 2 (1320)vr, / 2 (1270)<t, / 2 (1270)(7rvr)^ e , 

/ 2 (1270)p°, / 2 (1270)/ 2 (1270) 

7T*(1300)7T 



Three parameterisations of the a have been tried, that of Au, Morgan and Pennington [10 



that of Zou and Bugg |TTJ and the parameterisation found to fit the a contribution in the 



7r 7r°7r°7r 



channel described above. Three different parameterisations of the 7r*(1300) have been 
tried. In the Crystal Barrel analysis two parameterisations of the tt*(1300) have been found, 
one with M = 1114 MeV, T = 340 MeV M and a second with M = 1400 MeV, T = 275 MeV HI. 



The third parameterisation uses the PDG values |L2 



Different combinations of waves and isobars have been tried and insignificant contributions 
have been removed from the final fit. The best fit is shown in fig. |2| and consists of the following 
waves: J PC = 1++ pp with | J z \ = 1 fig. |b), J PC = 0++ pp fig. |c), J PC = 0++ a a fig. |d), 
jpc = 2 -+ a2 (i320)7r with | J z \ = 1 fig. |e) and J PC = 2++ f 2 (1270)nn with J z = fig. @f). 
As in the case of the analysis of the r]7i + 7i^ channel there is no need for any J PC = 2++ 
a2(1320)7r wave. The J PC = ++ 00 wave is only required if the parameterisation of the a 
found from the 7r 7T 7r 7T analysis is used. If the parameterisations of the a used to fit the 
7r + 7r~7r + 7r~ channel in our previous publication || are used here, the same conclusion would 
be drawn, i.e. no J PC = 0++ 00 wave is required. Hence as was stated in the introduction 
the parameterisation used to describe the a is crucial. There is no need for any J PC = 0++ 
7r*(1300)7r wave irrespective of the parameterisation used. 

Superimposed on the J p = 1 + pp wave shown in fig. ^|b) is a Breit Wigner convoluted with 
a Gaussian used to describe the /i(1285) in the fit to the 7r + 7r~7r + 7r~ mass spectrum [Q. As 
can be seen the /i(1285) is well described. 

The J p = + pp distribution in fig. 0c) shows a peak at 1.45 GeV together with a broad 
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enhancement around 2 GeV. A fit has been performed to the J p = + pp amplitude in fig. |2]c) 
using a single channel K matrix formalism |13|] including poles to describe the interference 



between the / (1370), the / (1500) and a possible state at 2 GeV. No account has been made 
for the pp threshold in this fit. The result of the fit is superimposed on the J p = + pp 
distribution shown in fig. |2]c) and well describes the data. The resulting T-matrix sheet II pole 



positions [14] for the resonances are 



/o(1370) M= (1309 ±24) -i (163 ± 26) MeV 
/o(1500) M= (1513 ±12) — i ( 58 ± 12) MeV 
/ (2000) M= (1989 ±22) -if 224 ±42) MeV 



These parameters are consistent with the PDG ]I2| values for the /o(1370) and /o(1500) 



The J p = + aa distribution in fig. ||d) shows a peak at 1.5 GeV. Superimposed on fig. |]d) 
is what would be expected from the / (1500); as can be seen there is little evidence for any 
/o(1370) contribution. A result similar to what was found in the analysis of the 7r 7r°7r 7r° final 
state. Correcting for the unseen pp and aa decay modes the branching ratio of / (1500) to 
pp/aa = 3.3 ± 0.5. 

These results disagree with the values found by the Crystal Barrel collaboration 0, || which 
found a considerably larger contribution of aa in the decays of the /o(1370) and /o(1500). In 
order to determine the maximum aa contribution the parameters of the a would need to be 
changed such that it has a similar shape to the p. In this case the difference between a pp 
and aa decay mode would be determined by the relative angular momentum between the pions 
from the decay of the p or a. Using such a parameterisation for the a we can determine a 
maximum contribution in the decay of the / (1370) to aa of 30 % of the pp decay. For the 
/o(1500) such a parameterisation of the a would yield equal decay rates to aa and pp. 

As in the case of the r]7i + ii' - channel 0, the J PC = 2~ + a 2 (1320) 7r wave is consistent with 
being due to two resonances, the 772 (1645) and the 772 (1870). Superimposed on fig. ^e) is the 
result of a fit using two resonances to describe the 772(1645) and 7/2(1870). The masses and 
widths determined for each resonance are given in table [I] and are consistent with those found 
previously in the r/im final state ||. 

Fig. |2|f) shows the J p = 2 + /2(1270)(7T7r) S ' wai , e distribution which has been fitted with a 
single Breit-Wigner, with M = 1980 ± 22 MeV and T = 520 ± 50 MeV. These values are 
compatible with those coming from the fit to ^(1950) from the 7r + 7r~7r + 7r _ mass spectrum Q. 

Finally, a reanalysis of the 7r + 7r~7r + 7r~ channel has now been performed using the new 
parameterisation of the a described above. A new version of the Monte Carlo program has also 
been used which has been improved since the previous publication Q. The 7r + 7r~7r + 7r~ effective 
mass spectrum is shown in fig. |3]a) and consists of 1 167 089 centrally produced events. 

The analysis has been performed independently from the 7r + 7r _ 7r 7r° analysis and different 
combinations of waves and isobars have been tried and insignificant contributions have been 
removed from the final fit. As in the original analysis [§J of the 7r + 7r~7r + 7r~ system the addition 
of the J PC = 2 ++ a2(1320)7r wave improved the Log Likelihood in the 1.8 to 1.9 GeV mass 
interval and has the effect of splitting the / 2 (1950) signal. If we use the information from the 
?77r + 7r~ || and 7r + 7r~7r°7r° analysis and reject this wave then the best fit is shown in fig. |3] and 
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consists of the following waves: J PC = 1 ++ pp with \J Z \ = 1 fig. |b), J PC = ++ pp fig. |]c), 
jpc = ++ aa fig> jpc = 2 -+ a 2 (i320)7r with | J z \ = 1 fig. |e) and J PC = 2++ / 2 (1270)7T7t 
with = fig. |3|f). The results are then very similar to those found for the 7r + 7r~7r°7r° channel. 

A fit to the J p = 2+ f 2 {l270){7T7T) Sw ave distribution shown in fig. |f) gives M = 1940 ± 22 
MeV and T = 485 ± 55 MeV. In order to get further confidence that the /2(1950) has a dominant 
/2(1270)7T7r decay mode a measurement of its rate to 7r + 7r~7r + 7r~ , 7r + 7r~7r°7r° and 7T°7r 7r 7r has 
been made. The predicted ratio for an I = state would be 7r + 7r~7r + 7r~ :7r + 7r~7r 7r° :7r°7r°7r 7r = 
4:4:1. The measured rate is 4 : 3.7 ±0.2 : 0.95±0.1 consistent with the prediction. However, 
if the a2(1320)7r wave was introduced in the 7r + 7r _ 7r + 7r~ analysis the measured 7r + 7r~7r + 7r~ rate 
would drop by a factor of 2. 

Therefore, although the Log Likelihood would indicate the presence of a J PC = 2 ++ a 2 (1320)7r 
wave in the fit of the tt + 7t~7t + 7t~ channel, the lack of it in the 7/7r + 7r~ and 7r + 7r~7r°7r° together 
with the fact that the measured decay rates would be wrong indicates how dangerous it can 
be to rely on Log Likelihood from a single channel alone to determine the existence of a decay 
mode. 

Superimposed on the J p = 1 + pp wave shown in fig. [5]b) is a Breit Wigner convoluted with 
a Gaussian used to describe the /i(1285) in the fit to the mass spectrum. As can be seen the 
/i(1285) is well described. 

A fit to the J p = + pp distribution shown in fig. |3]c) has been performed using the parametri- 
sation used in the 7r + 7r~7r°7r° channel. The T-matrix sheet II pole positions for the resonances 
are 

/o(1370) M= (1295 ±24) -z(169±26) MeV 
/o(1500) M = (1509 ± 12) -i ( 44 ± 12) MeV 
/o(2000) M= (1995 ±22) -i(218±42) MeV 

These parameters are in good agreement with those found in the 7t + tt~7t 7t final state. 

The J p = + aa distribution in fig. |3]d) shows a peak at 1.5 GeV. Superimposed on fig. |]d) 
is what would be expected from the /o(1500). Similar to what was found in the 7r°7r 7r 7r° and 
7r + 7r~7r°7r° final states there is little evidence for any /o(1370) contribution. Correcting for the 
unseen pp and aa decay modes the branching ratio of / (1500) to pp/uu = 2.6 ± 0.4 which is 
consistent with the value found from the 7r + 7r~7r°7r° final state. As in the 7r + 7r~7r°7r° channel 
a parameterisation of the a has been used with a shape similar to that of the p(770) which 
give an upper limit to the aa contribution in the decays of the /o(1370) and /o(1500). For the 
/o(1370) the maximum contribution is 25 % of the pp decay. For the /o(1500) the maximum 
aa contribution is 80 % of the pp decay. 

Fig. ^e) shows the J PC = 2~ + a 2 (1320)7r wave. Superimposed is the result of a fit using two 
resonances to describe the 7/2(1645) and 7/2(1870). The masses and widths determined for each 
resonance are given in table [TJ and as can be seen are consistent with those found for the rjnn M 
and 7r + 7r~7r°7T final states. Table [l| also gives the combined masses and widths from all the 
decay modes, including those found in the 7/7r + 7r~ channel 0, where the common systematic 
errors have been taken into account. In both the 7r + 7r~7r°7r° and 7r + 7r~7r + 7r - channels there is 
also evidence for a J p = 2 + pp wave. Except for evidence for the ^(1270) this wave is a broad, 
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structure-less distribution. 

In order to gain further confidence in the results, we have measured the rate to 7r + 7r~7r + 7r~ , 
7r + 7r _ 7r°7r° and 7r°7r 7r 7r° for the other waves observed. As can be seen from table ^| the results 
obtained are consistent with the expectations. 

In the analyses presented in this paper, the decay modes of the /o(1370) and /o(1500) are 
in disagreement with what has been found by the Crystal Barrel collaboration. In order to 
demonstrate the clear need for a sizable pp contribution to the peak at 1.45 GeV a study has 
been made of the 2ir and 37r effective mass spectra for the mass range 1.4 < M(47r) < 1.5 GeV 
and for <1Pt < 0.2 GeV for which the peak at 1.45 GeV is most prominent [§J. Fig. f|a) and b) 
show the ir + 7T~ and 7r + 7r _ 7r ± mass spectra respectively from this region for the 7r + 7r~7r + 7r~ final 
state. Superimposed on the mass spectra as a dashed line is what would be expected from a pp 
final state and shaded is the contribution from a aa final state using the parametrisation of the 
a found from the 7r 7r 7r°7r° channel. The solid curve shows the sum of the two contributions and 
as can be seen, apart from an excess of events in the K° mass region of the 2tt mass spectrum, 
the simulation well describes the experimental tt + tt" and 7r + 7r _ 7r ± mass spectra. Therefore it 
can be seen that a strong pp contribution is needed in this mass region and there is no need 
for any 7r*(1300) contribution to the 3tc mass spectrum. A similar conclusion is found from an 
analysis for the 7r + 7r - 7r°7r° final state. 

In summary, there is evidence for an a2(1320)7r decay mode of the r/2(1645) and 7/2(1870) 
in the 7r + 7r - 7r°7r° and 7r + 7r~7r + 7r~ final states. The /2(1950) is consistent with being a single 
resonance with a dominant /2(1270)7T7t decay mode. The /o(1370) is found to decay dominantly 
to pp while the /o(1500) is found to decay to pp and aa. 
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Tables 

Table 1: Parameters of the r? 2 (1645) and r? 2 (1870) 



Resonance 


Final state 


Mass (MeV) 


Width (MeV) 


?? 2 (1645) 
^(1645) 


7r + 7r~7r°7r° 


1621 ± 11 
1624 ± 12 


168 ± 23 
173 ± 22 


^72(1645) 


Combined 


1617 ± 8 


177 ± 18 


r?2(1870) 
^(1870) 


7r + 7r~7r°7r° 


1862 ± 19 
1843 ± 23 


222 ± 28 
226 ± 32 


^(1870) 


Combined 


1844 ± 13 


228 ± 23 
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Table 2: The predicted rate for I = and measured rate to it + it it + it : it + it it it : 
7r 7r°7r 7r° for different waves. 



Wave 


Predicted ratio 


Measured ratio 


2++ / 2 (1270)7T7T 


4 


4 


1 


4 : 3.7±0.2 : 0.95±0.2 


1++ pp 


1 


2 





1 : 2.1±0.2 : 0. 


0++ pp 


1 


2 





1 : 1.9±0.2 : 0. 


2~+ a 2 (1320)7r 


1 


2 





1 : 2.0±0.2 : 0. 


0++ aa 


4 


4 


1 


4 : 3.9±0.2 : 1.2±0.3 
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Figures 



Figure 1: The 7r°7r 7r 7r° channel, a) M(77) when the other three 77 pairs lie in the n° region, b) 
M(7r 7r°7r 7r° ) and c) M(n°7r°) for 1.3 < M(7r 7r 7r 7r° ) < 1.5 GeV. Results of the spin analysis 
d) the J PC = ++ a a wave and e) the J PC = 2 ++ / 2 (1270)7T7r wave with fits described in the 
text. 



Figure 2: The 7r + 7r~7r°7r° channel, a) The total mass spectrum, b) 1 ++ pp, c) ++ pp, d) ++ aa, 
e) 2~ + a2(1320)7r and f) 2 ++ / 2 (1270)7r7r. The superimposed curves are the resonance contri- 
butions coming from the fits described in the text. 



Figure 3: The 7r + 7r~7r + 7r~ channel, a) The total mass spectrum, b) 1 ++ pp, c) ++ pp, 
d) ++ (jo-, e) 2~ + a 2 (1320)7r and f) 2 ++ / 2 (1270)7T7r. The superimposed curves are the resonance 
contributions coming from the fits described in the text. 



Figure 4: a) The 7r + 7r~ and b) ii + 7i^7i ± mass spectra for the 7r + 7r~7r + 7r~ final state in the 
mass range 1.4 < M(7r + 7r~7r + 7r~ ) < 1.5 GeV and with cIPt < 0.2 GeV. The superimposed 
histograms are for a pp (dashed) and era (shaded) final state. The solid line represents the sum 
of the pp (90 %) and aa (10 %) contributions. 
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